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(Dated: February 6, 2008) 

By studying the structures of clusters bound by a model potential that favours polytetrahedral 
order, we find a previously unknown series of 'magic numbers' (i.e. sizes of special stability) 
whose polytetrahedral structures are characterized by disclination networks that are analogous 
to hydrocarbons. 

PACS numbers: 61.46.+w,36.40.Mr 



Polytetrahedral order fjL |) has become an increasingly 
important concept in condensed-matter physics. Such 
polytetrahedral structures, for which the whole of space 
can be naturally divided up into tetrahedra with atoms 
at their vertices, are the basis of crystalline Frank-Kasper 
phases || Q, and have been invoked in order to under- 
stand the structure of quasicrystals [|[ ||, Q and atomic 
liquids and glasses ||, However, little is known 
about the consequences of polytetrahedral order for the 
structure of clusters and nanoparticles. This situation 
contrasts with close-packing schemes, which give rise to 
fascinating cluster structures, such as Mackay icosahedra 
fiof , Marks decahedra jllj and Leary tetrahedra 0|. 

However, recent experiments indicate that small cobalt 
clusters can have polytetrahedral order jl3j, although 
only limited information about their detailed structure 
could be obtained. Furthermore, there is an increasing 
interest in mixed metal clusters Jl4|, |l5), and polytetra- 
hedral structures would be expected for those alloys that 
exhibit Frank-Kasper or quasicrystalline phases in bulk. 

The distinctive features of polytetrahedral packings 
stem from the inability of regular tetrahedra to fill all 
space. When five regular tetrahedra are packed around 
a common edge there is a small angular deficit of 7.4°, 
whose closure requires a small distortion of the tetrahe- 
dra. If this method of packing is extended to larger collec- 
tions of tetrahedra, local icosahedral coordination results, 
but the strain that needs to be introduced to close all the 
gaps grows very rapidly. Therefore, in order to form an 
extended polytetrahedral structure, sites where six tetra- 
hedra share a common edge need to be introduced — a 
negative disclination line is said to run along this com- 
mon edge. Even though the local distortion required to 
remove the overlap that occurs when packing six regualar 
tetrahedra is larger, the overall strain is reduced. 

Polytetrahedral packings can therefore be described by 
a network of disclination lines threading an icosahedrally- 
coordinated medium. In Frank-Kasper crystals this 
disclination network is ordered and periodic, whereas it 
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FIG. 1: Comparison of the Dzugutov potential in its orig- 
inal (Dz) and modified (m-Dz) form. 



has been suggested that atomic liquids and glasses are 
characterized by disordered entangled disclination net- 
works n §. 

Many very small clusters naturally form polytetra- 
hedral clusters. The 13-atom icosahedron (20 slightly 
distorted tetrahedra sharing a vertex) is an extremely 
common structure for rare gases metal |0 and 

molecular Jl8| clusters and is generally favoured over a 
close-packed structure because of its lower surface en- 
ergy. However, structures that continue this polytetrahe- 
dral packing soon become disfavoured for most materials 
because the associated strain cannot be accommodated. 
The largest polytetrahedral clusters have been obtained 
for a model system where the width of the potential al- 
lows the system to tolerate ordered polytetrahedral struc- 
tures up to 7V«70 

A model system that exhibits polytetrahedral clusters 
and enables us to make structural predictions about such 
clusters would therefore be of great interest. Here we seek 
to address this issue by studying clusters interacting with 
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a potential of the Dzugutov form 1 20 



V(r) = A(r- m - B)exp 



Q(a -r) + 



Bexp 



r — b 



6(6 -r), 



(1) 



where 0(x), the Heaviside step function, is for x < 
and 1 otherwise. The total potential energy of a clus- 
ter is then E — J2i<j V( r ij)> where is the distance 
between atoms i and j. This potential was designed to 
encourage polytetrahedral and local icosahedral order in 
supercooled liquids through the introduction of a local 
maximum in the potential near to \[2 times the equilib- 
rium pair distance (Fig. that disfavours close-packed 
structures |plf . This maximum also somewhat resem- 
bles the first of the Friedel oscillations that can occur for 
metal potentials. Interestingly, under certain conditions 
this potential is able to form a dodecagonal quasicrystal 
II- 

For the original parameterization of the Dzugutov po- 
tential, clusters form non-compact polytetrahedral struc- 
tures composed of needles, rings and three-dimensional 
networks of interpenetrating and face-sharing icosahcdra 
p3| . Thus the original potential cannot provide a realis- 
tic model of the compact polytetrahedral clusters formed 
for cobalt or that might occur for metallic alloys. Non- 
compact structures occur because the relatively narrow 
potential well (Fig. [l]) does not allow the system to ac- 
commodate the strain associated with compact polyte- 
trahedral clusters. Therefore, we chose a new parame- 
terization of the Dzugutov potential that both increases 
the width of the minimum and reduces the width of the 
potential maximum. As envisaged, this new potential 
gives rise to compact polytetrahedral clusters. The pa- 
rameters in our modified potential have the values 



A = 3.00 5 = 2.109 a = 1.65 6=1.94 
c = 0.52 d = 0.55 m = 4. 



(2) 



The pair potential has a maximum at r max = 1.36 r cq of 
height 0.83e, where r eq and e are the equilibrium pair 
separation and well depth, respectively. 

To find the global minima of clusters interacting 
with this potential we used the basin-hopping algorithm 
p4] , p5[ . For each size up to 7V=100 we performed five 
runs of 100 000 steps starting from a random configura- 
tion. We also performed short runs starting from con- 
figurations generated by adding or removing the appro- 
priate number of atoms from some of the lowest-energy 
minima for sizes one, two or three atoms above or be- 
low the current size. These latter seeded runs were re- 
peated until no new putative global minima were found 
and were particularly important because the roughness of 
the energy landscape |2^| makes optimization from a ran- 
dom starting point particularly difficult. Indeed above 
7V=100 no runs from a random starting point found the 
global minima, so in this size range we used a differ- 
ent approach. Using the structural principles obtained 
from the particularly stable structures for clusters with 
less than 100 atoms, we were able to construct a series 
of candidate geometries for particularly stable sizes in 
the range iV=100-250. These structures then served as 
the initial seed configurations for a series of short basin- 
hopping runs for the intervening sizes. Again seeded runs 
using the lowest-energy minima of nearby sizes were re- 
peated until no further improvements were obtained for 
any size. 

The energies of the resulting putative global minima 
are depicted in Fig. |^ in a manner that emphasises par- 
ticularly stable minima or 'magic numbers'. At small 
sizes the magic numbers (AT=13, 19, 23, 26, 29) are those 
expected for polytetrahedral growth upon the 13-atom 
icosahedron and are similar to those seen for Lennard- 
Jones clusters 
gon clusters \U7 



26J and their experimental analogue, ar- 
The next magic number corresponds to 
the disc-like 38-atom cluster, which was previously found 
for clusters interacting with the original Dzugutov poten- 
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FIG. 3: Structures of some of the magic number clusters. 
On the right is the complete structure, in the middle the 
corresponding disclination network and on the left the 
structure that is at the centre of the cluster. All three 
have the same orientation. 



tial and which has a single disclination line running along 
its axis |23| . Then, for the rest of the size range we con- 
sider, there are series of roughly equally spaced minima 
in Fig. ^| that correspond to a new sequence of magic 
numbers. 

Some of the structures of these new clusters are de- 
picted in the right-hand column of Fig. ||. Each clus- 
ter consists of a disclinated central structure that is sur- 
rounded by an overlayer in which an atom is added to 
each face and above each vertex that is not at the end of a 
disclination line. This is the same overlayer as for the ini- 
tial polytetrahedral growth on the 13-atom icosahedron, 
and in that context has been called the anti-Mackay |n]] 
or face-capping |26j| overlayer. This overlayer does not ex- 
tend the disclinations of the central cluster, thus giving 
rise to characteristic six-fold pits where the disclinations 
exit onto the surface of the resultant cluster (Fig. ||). 

At the centre of the 57-atom cluster is a 17-atom struc- 
ture in which the central atom is the node for four discli- 
nations in a tetrahedral arrangement. This coordination 
shell is commonly found in the Frank-Kasper phases [Q. 
Furthermore, the atoms of all the central clusters either 
have this or an icosahedral coordination shell, thus giving 
rise to a tetrahedrally coordinated disclination network 



(Fig. |J). There are two bulk Frank-Kasper phases that 
involve such networks, the C14 and C15 phases. In the 
C15 phase the disclination network has the structure of 
the diamond lattice and in the C14 phase the wurtzite 
structure. All the cluster centres can be considered to be 
fragments of these two Frank-Kasper phases, except the 
centre of the 221-atom cluster which involves a mixture 
of the two phases. 

The easiest way to understand the progression of struc- 
tures is to note the correspondence between the discli- 
nation networks and hydrocarbon structures (Fig. ||). 
The disclination network of the 57-atom cluster is anal- 
ogous to methane. Next comes a series corresponding 
to the linear alkanes, ethane (76), propane (94) and n- 
butane (112), and the branched alkane, isobutane (111). 
At this point it becomes favourable to form more com- 
pact structures analogous to cycloalkanes, e.g. the chair 
form of cyclohcxanc (129) and methyl-cyclohexane (146). 
Above this size structures analogous to cage hydrocar- 
bons Q are favourable, such as bicyclo[2.2.2]octane 
(148), adamantane (166), diamantane (203) and triaman- 
tane (239). It is noticeable that the most stable of these 
latter magic numbers correspond to the polymantanes 
(or diamondoids), where the central structures are frag- 
ments of the C15 phase. C14 and mixed disclination net- 
works are only competitive in between these sizes when 
they are competing with structures analogous to methyl- 
polymantanes. 

The analogy to the hydrocarbons also allows ready pre- 
diction of the structure of larger clusters. For example it 
is well-known that the next polymantane, tetramantanc 
has three isomers j2£|. However, as for the original pa- 
rameterization of the Dzugutov potential |2^] , the body- 
centred-cubic (bcc) lattice is lowest in energy for bulk. A 
comparison of E ave to a similar function fitted to a series 
of bcc rhombic dodecahedra indicates that bcc clusters 
are lowest in energy beyond ArT400. 

Experiments can usually only probe cluster structure 
indirectly and so comparisons with calculations from can- 
didate geometries are required for structural identifica- 
tion. Therefore, to enable the potential identification of 
the structures described here, all are accessible from the 
Cambridge Cluster Database (30). Furthermore, as the 
series of magic numbers does not coincide with any pre- 
viously known fT7j| , there is also the potential for identi- 
fication through mass spectral abundances. 

Calculations for an example cluster, the 166-atom 
adamantane analogue, are consistent with the experimen- 
tal results for cobalt clusters fL3[ . The calculated scat- 
tering function reproduces many of the features observed 
in the x-ray data, and the structure has perpendicular 
planes of atoms that can account for the square arrays 
of lattice fringes seen in the electron microscopy images. 
However, this agreement may result more from the poly- 
tetrahedral order than the actual detailed structure of 
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FIG. 4: Schematic representations of the disclination net- 
works for the particularly stable clusters. The terminal 
disclination lines have been omitted for clarity. The net- 
works, thus drawn, are analogous to the carbon back- 
bones of a series of hydrocarbons. 



our clusters |p.3| . 

It is remarkable that many of the unusual cluster struc- 
tures that have been predicted by theory have been sub- 
sequently observed experimentally. For example, the 
Mackay icosahedron (To), first suggested in 1962, has 
since been seen for a wide variety of systems over a 
large size range [jlT], |3l]| . Furthermore, the small trun- 
cated octahedron at N—38 and Marks decahedron at 
iV=75, whose stability was first identified in calcula- 
tions on model clusters |32|, have since been identified 
for Ni 38 Q, Au 38 @], and Au 75 ||. Most recently 
the Leary tetrahedron, a surprise global minimum for a 
98-atom Lennard- Jones cluster 12 , has since been found 
for (C6o)98 [^6|- Therefore, it would be no surprise if the 
polytetrahedral structures that we have described here 
were likewise to be positively identified experimentally. 

JPKD is grateful to Emmanuel College, Cambridge for 
financial support. 
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